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Abstract 

We investigate the percolation properties of a two-state (occupied - empty) 
cellular automaton, where at each time step a cluster of occupied sites is 
removed and the same number of randomly chosen empty sites are occupied 
again. We find a finite region of critical behavior, formation of synchronized 
stripes, additional phase transitions, as well as violation of the usual finite- 
size scaling and hyperscaling relations, phenomena that are very different 
from conventional percolation systems. We explain the mechanisms behind 

all these phenomena using computer simulations and analytic arguments. 
PACS numbers: 64.60.Lx, 05.70.Jk, 05.70.Ln 
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I. INTRODUCTION 



During the past years, systems which exhibit self-organized criticality (SOC) have at- 
tracted much attention, since they might explain part of the abundance of fractal structures 



sand grains M, increase of strain 0, tree growth ||, spontaneous mutations M) and rare 
dissipation events which are instantaneous on the time scale of driving (e.g. sand avalanches, 
earthquakes, fires, or a series of rapid mutations). In the stationary state, the size distribu- 
tion of dissipation events obeys a power law, irrespective of initial conditions and without 
the need to fine-tune parameters. There is, however, no reason to expect that systems with 
slow driving and instantaneous avalanches are always SOC. Such systems might also have 
many small avalanches which release only little energy, or only large avalanches which re- 
lease a finite part of the system's energy, or some combination of both. SOC systems are 
naturally at the critical point, due to e.g. a conservation law (sandpile model), a second 
time scale separation (forest-fire model), a competition between open boundary conditions 
and the tendency of neighboring sites to synchronize (earthquake model HQ - see however 
H for a counterexample), or due to extremal dynamics ('evolution' model Q). Often, the 
critical behavior breaks down when details of the model rules are changed (e.g. the bound- 
ary conditions in the earthquake model || or the tree growth rule in the forest-fire model 



There are certain parallels between these models and equilibrium critical systems, since 
both consist of many small units which interact with their neighbors, and since spin clusters 
in an Ising model or clusters of occupied sites in percolation theory can be compared to 
avalanches. However, the critical behavior of nonequilibrium systems can depend on mi- 
croscopic details, as mentioned above, in contrast to equilibrium critical phenomena, which 
commonly show universal behavior. Also, nonequilibrium systems do not satisfy a detailed- 
balance condition and can e.g. show periodic behavior. Furthermore, avalanches are usually 
released when some variable reaches locally a threshold, while other regions of the system 



in nature M . Their common features 
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might be far below the threshold, and consequently not all parts of the system look equal. 
This can in particular result in more than one diverging length scale, as in the earthquake 
model H or in the forest-fire model ||. By contrast, in equilibrium systems the energy is 
an extensive variable, which means that all regions (which are large compared to the lattice 
constant) are equal. This is e.g. the basis for hyperscaling relations. 

In this paper, we discuss in detail a model |§ that belongs to the mentioned class of 
nonequilibrium systems with avalanche-like dynamics. It is a nonequilibrium percolation 
model, where clusters of occupied sites are removed, and the same number of sites that 
have become empty are occupied at random. The density of occupied sites is the control 
parameter of the model. The 'avalanches' of our model are removal events, and the size of 
an avalanche is the size of a removed cluster. This model illustrates well the fundamental 
differences between equilibrium and nonequilibrium, showing various features that are not 
observed in equilibrium systems: The region of small avalanches and the region of infinite 
avalanches are separated by a finite region of critical behavior, where the correlation length 
diverges slower than the system size. The exponent that relates the system size with the 
correlation length depends on the density. Besides the correlation length, there are other 
relevant length scales. Since the critical behavior occurs over a finite density interval, the 
system can exhibit power laws naturally, without fine tuning of parameters to a precise 
value. Therefore, our model belongs to the class of SOC systems. In the region of infinite 
avalanches, the system shows synchronization with a period that depends on the value of 
the density. We illustrate and explain all these observations using computer simulations and 
analytical arguments. Part of the results were already published in 0. 

The work is structured as follows: In Sec. J|, we define the model. In Sec. [II, the 



subcritical phase and the critical point of the model are treated. The mechanism that leads 
to criticality and the value of the critical density are explained, and the exponent of the 
cluster size distribution in ID is calculated analytically. Sec. [TV] discusses the critical phase. 
The reason for the existence of a whole critical phase as well as for its properties like non- 
standard finite-size scaling and violation of hyperscaling are explained. The supercritical 



phase is treated in Sec. [V]. First, we explain the existence of synchronized stripes and their 
relation to the subcritical phase (Subs. |V A| ), then we discuss hysteresis and the maximum 
possible number of stripes (Subs. |V ij| ) , and, finally, we investigate stability, movement and 
roughness of stripes (Subs. | V C| ) . In the conclusion Sec. [VT], we summarize and discuss our 
work. 

II. DEFINITION OF THE MODEL 

The non-equilibrium percolation model model is defined on a <i-dimensional hyper-cubic 
lattice with L d sites. Each site is either occupied or empty. The control parameter of the 
system is the density of occupied sites, p. 

The dynamics are defined by the following rules: (i) An occupied site in the system is 
chosen at random and the whole cluster of s occupied sites connected to this site (by nearest- 
neighbor coupling) is removed from the system, i.e. the occupied sites of that cluster turn 
to empty sites, (ii) We occupy s randomly chosen empty sites (possibly also including sites 
which have become empty due to the removal of the cluster), (iii) Proceed with (i). 

These rules ensure that the density of occupied sites p is a conserved quantity. Starting 
with a random initial state, the system approaches after a transient time a stationary state 
that is characterized by a certain size distribution of clusters, and where time average and 
ensemble average of all quantities are identical. Throughout this paper, we discuss only the 
properties of the stationary state. These properties, that are explained in detail in the next 
three sections, are as follows: For small densities, there are only small clusters of occupied 
sites in the system. With increasing density, the size of the largest cluster increases, and it 
diverges at a critical density p c . For p c < p < p?' , the system is critical, i.e. the cluster size 
distribution is a power law. The size of the largest cluster diverges slower than the system 
size. For p > p^\ the system has a finite number of regions of different density. The region 
with the highest density has a spanning cluster. 

One can think of the dynamics of this model as 'explosions' that take place at a randomly 
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chosen site. During the explosion, the whole cluster connected to the explosion site is blown 
up, and its constituents settle down somewhere else in the system. Or one might think of 
colonies of animals which are dispersed into all directions by some enemy or other event. 
From a more abstract point of view, one has a nonequilibrium percolation problem. 

This model is also closely related to the self-organized critical forest-fire model (SOC 
FFM) ID, when occupied sites are equated with trees: For p < p c , the correlation length 
£ and the mean number of removed sites per step, S, are finite, and all properties of the 
stationary state can be found by looking at a section of the system of linear size £, with 
( < f C I. In this section, there is no conservation of density, and the dynamics can 
be characterized by a small 'tree growth' rate p = Sp(£/L) d /(1 — p), and a 'lightning' 
rate / = (£/L) d . The tree growth rate is the probability that a given empty site becomes 
occupied during one step, and the lightning rate is the probability that a given site is 'struck 
by lightning' (i.e., selected) per step, with the consequence that all 'trees' connected to this 
site 'burn down' (i.e., are removed). S diverges in the limit f/p—> as 

S = p(l-p)/fp. (1) 

In this limit, p approaches the critical value p c . The dynamics in the small region of size 
I are the same as for the SOC FFM, and the critical exponents close to p c are therefore 
the same as those of the SOC FFM in the limit f/p — > 0, when the //p-dependence is 
translated properly into a p-dependence (simulation results of the SOC FFM can be found 
in e.g. [p!0|-|T3|,|8[|). Choosing p as the control parameter instead of f/p, allows us to study 
the FFM beyond the critical point. This was in fact our original motivation to introduce 
the model studied in this paper. 

In the following three sections, we discuss in detail the subcritical, the critical, and the 
supercritical behavior of the model. Unless stated otherwise, the considered system is a 
two-dimensional square lattice. 
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III. SUBCRITICAL PHASE AND APPROACH TO THE CRITICAL POINT 



First, we discuss the parameter region p < p c , where the system has a cutoff in cluster 
size that is independent of the system size L. Let S be the mean number of sites removed 
from the system in one time step, without taking into account the refilling of sites. For very 
small densities, there exist only very small clusters, and most clusters will consist of only 
one site. The process of removing clusters (in this case mostly single sites) and refilling sites 
at random does not change this situation, i.e. S is of the order of one, and the occupied 
sites remain randomly distributed, as for a percolation system at small densities. 

With increasing density, S increases also. Removing the chosen cluster and refilling its 
sites at random into the system introduces fluctuations in the local density of occupied sites, 
because the removed cluster leaves behind a 'hole' and the refilled sites increase the density 
in the rest of the system to a value larger than p. The critical density p c , where S diverges, 
and the critical exponents close to p c , are therefore different from their values in percolation 
theory. 

We define the usual quantities that are investigated in percolation systems (for an intro- 



duction to percolation, see [[TJJ): The number density of clusters of occupied sites of size s 
will be denoted by n(s). Near a critical density p c , n(s) is expected to behave like a power 
law 

n(s) oc s~ r C(s/,s max ), (2) 

where C is a scaling function and s max oc \p c — p|~ 1 / CT . The average cluster size S is defined 
by 



and is expected to diverge like 



S oc \p c - p\ 7 . (4) 
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The correlation length £ is defined as the root mean square distance between occupied sites 
on the same cluster, averaged over all clusters, which leads to 



. 2 ET=iR 2 (s)s 2 n(s) 



Near a critical point, £ is expected to diverge like 

£oc| Pc - P r. (6) 

The radius of gyration i? of a cluster grows with its size s like 

R(s) oc (7) 

with the fractal dimension df. These critical exponents are related via the scaling relations 
1/(7 = 7/(3 — r) = dfV. Finally, the strength of an infinite cluster is denoted by P. In 
percolation, P follows a power law 

P OC (p - Pc f 

above the critical point. 

In our simulations, we found p c ~ 40.8%, r = 2.15(3), df = 1.96(2), v = 1.20(5), and 
7 = 2.09(5). The values of these exponents as well as the value of the critical density p c are 
different from percolation theory, and they are identical with the corresponding values of the 
SOC FFM. The (p c — p)-dependence of our model can be translated into the //p-dependence 
of the FFM using Eq. (|j) and Eq. ([[]), giving 

p c -pcx U/Vf h - 

Thus, our values of a and v can be calculated from those of the FFM by multiplication 
with 7. r and df are exponents related to the cluster size s, so no multiplication with 7 is 
necessary. The exponent (3 vanishes, as we shall see in the next section. 

The above-mentioned fluctuations in the local density of occupied sites can easily be 
seen by looking at a snapshot of the system for densities close enough to the critical density. 



A typical stationary state for the density p = 39.3% and system size L = 1024 is shown 
in Fig. [I]. One can see that the system consists of a large number of regions with different 
and rather homogeneous density. The typical size of these 'patches' does not depend on 
the system size L, provided that L is large enough. Many properties of the model can 
be understood by describing the system in terms of these patches of homogeneous density 
of occupied sites. For small average patch size (like in Fig. |l|), it is not always possible 
to assign a given site unequivocally to a certain patch. This changes, however, when the 
critical density is approached, where the mean patch size is larger and the patch boundaries 
become sharper. 

For large system size L 2 ^> S, only a few sites are occupied in a given patch per time 
step, and the density p{t) in a patch evolves continuously according to p'(t) = p(l — p(t)) 
with some growth rate p, which leads to 

pit) =i_(i_ exp(-pt) . (8) 

Here, time is measured since the removal event that produced that patch, leaving behind 
a small density of occupied sites p(0) = p after . We also define the mean density p bcfore of a 
patch just before its spanning cluster is removed, and the time T that it takes to increase 
the density from p aftcr to p bcfore . We easily derive T = (l/p) ln((l — p aftcr )/(l - p before )). The 
average density of a patch is 

p = i f T dt(l - (1 - p^)e-^) 
1 Jo 

^before „aftcr 

= l _ L t (9) 

L1L \ l_pbcforc J 

We measured the average values p before 62.5% and p after 7.8%, similar to the values 
found in [|T5| ,|S|. Interestingly, the same values will play an important role in the critical phase 
and in the striped phase discussed in the following two sections. With the measured values 
of p before and p after ; we obtain as average density of one region p = 39.2%. For large enough 
system size L 2 and neglecting the interactions between the different regions, the time average 
of the overall density p would also be p = 39.2%. In a real system with interacting regions, 
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the mean density is in general different due to the following two mechanisms: (i) Temporal 
oscillations in patch size: During the growth process of the density of a patch from p after to 
p he , all the neighboring patches have one removal process on an average. During each of 
those removals, also some sites of the central patch are removed, because patch boundaries 
are not cluster boundaries. This leads to a shrinking of the area of the original patch while 
its density increases, reducing the mean density from the above calculated value. (Of course, 
the original patch size is ultimately restored when the spanning cluster of the central patch, 
and with it some finite clusters of neighboring patches, are finally removed.) The relative 
shrinking of the area is stronger for smaller mean size of removed patches, since the number 
of sites removed from the central patch depends on the length of its boundary. For large 
mean size of removed patches, the ratio of initial patch size to final patch size approaches 
1, i.e. p = 39.2%, which explains why the density is smaller than the critical density when 
the mean size of removed patches (of the order S) is finite. (Fluctuations in the size of 
removed patches are complicated to consider, and do not affect our main conclusions.) (ii) 
The critical density itself is larger than the above value p = 39.2%, because the density 
at which a removal event takes place fluctuates around its mean value p before . Assuming a 
homogeneous distribution of occupied sites within a region, we measured the mean density 
p for various densities p before between 59.3% and 70%. The result is shown in Fig. ^. One 
can see that this function has a minimum for the measured average value of p before « 62.5%. 
Any fluctuation will therefore lead to an increase in the mean value of p. The fact that the 
density is so close to its minimum indicates a tendency of the system to maximize energy 
dissipation. To obtain the order of magnitude of the shift in p c due to fluctuations, we 
calculated numerically the mean density in the simulated interval [p c ;0.70] of p before . The 
result is 40.5(3)%, which is close to the correct critical density p c ~ 40.8%, indicating that 
the fluctuations in p before can indeed induce the observed shift in the density. 

These considerations show that the critical state of our model can to good approximation 
be interpreted as a combination of percolation systems of different densities, as also sug- 
gested in ||. The cluster size distribution n(s) is therefore the superposition of the cluster 
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size distributions of all the patches with their different densities between p after and p before ; 
distributed according to Eq. (§). 

Let us first consider the one-dimensional system. There, the percolation threshold is one, 
and even in dense regions all clusters are finite percolation clusters. Removal of a cluster 
leaves behind a string of empty sites. The system is therefore composed of strings of size 
S of different densities that represent different stages in the growth process of an empty 
string until it is completely filled. A string of density p contains a cluster of size s with 
probability p s (l — p) 2 . If the system size L is large enough, the time average equals the 
ensemble average, and the cluster size distribution should be given by 



n(s) oc / dt p(t) s (l-p(t)Y 
dt (l-e- pt ) s (e- pt ) 2 



o 



o 



oc / dt VI" \(-l) k e- kt e- 2t 



k=0 



kj y ' 2 + k 

1 



(s + l)(s + 2) 
s~ 2 for large s 



in agreement with the exact calculation of ||16fl . Similarly, the size distribution of hole clusters 
is found to be 



oo 



h( s ) oc / dt p(ty(i-p(t)y 
i 







'OC 



s(s + l)(s + 2) 
s~ 3 for large s , 



again in agreement with exact results [ 17| , [T6 1 



For dimensions higher than one, unfortunately, these calculations can not be carried out. 
There, the number of neighbors t of a cluster of size s does not only depend on s, but also 
on the shape of the cluster, and the number of different clusters with a given size s and 
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a given perimeter t is not known exactly. Furthermore, the integration over t does not go 
from zero to infinity, since the lowest and highest densities are not and 1, respectively. 
We can, however, determine numerically the cluster size distribution of a combination of 
percolation systems of different densities. We measured n(s) in 20 homogeneous systems 
with L = 2048 and different densities between p after and p before ; distributed according to 
Eq. (H). The total n(s) of all systems oc J2i=i n i( s ) is plotted in Fig. ||, together with the 
cluster size distribution of a percolation system at p = p c ~ 0.593. The exponent r ps 2.15 for 
a combination of percolation systems is indistinguishable from the cluster size distribution 
exponent of our model, showing that the latter is indeed generated by a superposition of 
percolation systems. Fig. [| does not show the bump observed in our simulations for large 
s. For large cluster sizes, the size distribution of patches with density above the percolation 
threshold affects the cluster size distribution. This effect is not present in Fig. [3], where there 
are no different patch sizes. It remains an open question whether the exponent r remains 
the same for much larger values of s than those accessible to simulations. At very large 
scales, the emerging dynamics of patches might become an important factor in determining 
the cluster size distribution. 

Another conclusion that can be drawn, is that r has to be always larger than or equal 
to the t of the corresponding percolation problem. n(s) results from a superposition of 
many percolation systems, so that apart from the critical percolation system, there exist 
also many systems with finite cutoff in the cluster size distribution. These systems increase 
the weight of the small clusters, and the slope of n(s) decreases, resulting in a larger r. For 
all dimensions and lattice types investigated so far (see for example |i~3|), this was true. 

IV. THE CRITICAL PHASE, HYPERSCALING, AND FINITE-SIZE SCALING 

For p > p c , one might expect the appearance of an infinite cluster which spans the whole 
system, as in percolation theory. However, the values of the critical exponents for p < p c 
already indicate that the situation in our system is very different from percolation. In con- 
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trast to percolation, the hyperscaling relation d = df{r — 1) is violated, which means that 
not every part of the system contains a spanning cluster at criticality [ID . This can also be 



seen easily by looking at the snapshot of a system at p ~ p c in Fig. In addition to large 
patches with high density, there exist also large patches with lower densities. Most of the 
patches have densities below the percolation threshold and contain many finite clusters. The 
patches with high densities (between p c and m 62.5%) contain not only finite, but also 'infi- 
nite' (i.e., spanning) clusters. In contrast to ordinary percolation, there is no homogeneously 
distributed set of large clusters that could join at p = p c to form the infinite cluster that 
spans the whole system. Rather, besides the largest cluster, the system contains many other 
large clusters that represent different growth stages of the largest cluster itself. The critical 
behavior of our model occurs over a finite interval [p£^(= Pc ~ 40.8%), p^(~ 43.5%)], where 
the cutoff in cluster size s max diverges, but slower than L 2 . The correlation length diverges 
slower than L. Since it is not possible to define a truly infinite cluster like in percolation in 
this phase, all clusters in the system contribute to the defining equations of s max , £ and S. 
We find 

Smax ocL^,eocL^,5ocL^ 

with p-dependent exponents 0i,2,3, while r and df remain unchanged. Tab. | shows the 
values of the exponents for different densities. Fig. |5] shows a snapshot of the system for 
p = 0.43. 

Fig. H] shows the size distribution of clusters sn(s) for different system sizes at fixed p 
(a) before and (b) after rescaling. With Eqs. (@) - ([/]) one can derive the scaling relations 
<pi = df(p2 and 03 = (3 — r)0i, which are well confirmed by the simulations. For p — > p^\ 
0i.2,3 approach the values df, 1, and 1/z/, respectively. 

The latter values are those that one would already have expected at p = p^p from finite- 
size scaling theory. In a critical system, a quantity X that scales as \T — T c \~ x cx ^ x ^ u in an 
infinite system, in a finite system is expected to obey the scaling form 
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with 

{const for x C 1 
x -xl" f or x > i 

such that for £ <C L, X oc |T — T c |~ x like in the infinite system, and for L <C £, X oc L*/ 1 '. 
The standard finite-size scaling exponents of s max , £ and 5 would therefore be 1/va = df, 
vjv — 1, and 1/za 

The main difference between our system and percolation concerning finite-size scaling is 
that in percolation a system of length Li arge without finite-size effects can be generated by 
putting together systems with finite-size effects of length L sma ii- In our model, a reorganiza- 
tion takes place when smaller systems are put together, since the smaller systems now can 
have fluctuations in their number of occupied sites that they could not have when they were 
isolated. Due to the inhomogeneous nature of the stationary state, these fluctuations are 
very large. 

We obtain a deeper understanding of the critical behavior of our model when we describe 
the system again in terms of patches of different densities. In the previous section, we found 
that the density within a patch grows continuously for p < p c and for sufficiently large 
system size. For not sufficiently large system size, and for p > p c , however, the finite system 
size leads to discontinuous jumps in the density within a patch. The non-continuous version 
of Eq. (D is (see Eq. (|IJ) below) 

^before ^after 

9 ~ 1 ~ 7l(((l - p before )/(l - pafter))-l/n - 1) • 

for the density of a system that takes n time steps to go from p after to p before . In Fig. [7], 
p is plotted as function of the number of patches n. One can see that a decrease in n 
raises the mean density. (Of course, this conclusion can also be obtained from an analytical 
calculation.) This effect supplements the two other mechanisms that affect the mean density 
presented in the previous section. For fixed p and changing L, these three mechanisms have 
to be in balance with each other. While the step size decreases with increasing L, tending 
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to decrease the density, the size of the large patches and the fluctuations in p before increase, 
cancelling the effect of the smaller step size on the density. 

So far, we have not yet discussed the possibility of fusion and splitting of patches. As 
soon as two neighboring patches have a density above the percolation threshold, they fuse 
and will be removed together. This effect must be balanced by a mechanism that splits 
patches. As long as the density of a patch is below the percolation threshold, neighboring 
removals move the boundary of this patch inwards, and a splitting into two patches occurs 
when opposite boundaries meet. A splitting can also occur when a finite cluster connecting 
two opposite edges of the patch is removed before the density in the patch reaches the 
percolation threshold. Since even large patches must split at the same rate at which they 
fuse with neighboring patches, the shape of patches cannot be round, but must be 'fingered', 
with necks of a width that does not depend on the patch size. The snapshots Figs. ^ and |5| 
show the fingered structure. It implies some characteristic length scale, the 'finger thickness', 
besides the lattice constant and the correlation length. This scale is a function of the density 
and the system size. The existence of several length scales was also observed in ||. 

V. SUPERCRITICAL PHASE 
A. Synchronized States 

At p = p^\ the correlation length finally becomes proportional to the system size, and 
the system has an 'infinite' cluster that contains a finite percentage (independent of L) of 
all occupied sites in the system. 

Fig. |] shows a snapshot of the system for p = 0.45. We see a system with five homoge- 
neous and equally large stripes with different densities. The stripe with highest density is 
above the percolation threshold and contains an infinite cluster (with non-zero strength P) 
as well as some small clusters. The other stripes are below the percolation threshold and, 
consequently, contain only small clusters. When one of the small clusters in this state is 
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removed, only few sites are redistributed, and the state of the system remains essentially 
unchanged. When the infinite cluster is removed, a large portion of all occupied sites in the 
stripe with the highest density are redistributed all over the system. The stripe which used 
to have the highest density now has the lowest density, while the density of the other stripes 
has increased. The values of the five densities are the same as before, except that they are 
now associated with different stripes. If we measure time in units of large redistributions of 
sites, the state of the system is periodic with period 5. 

Increasing p, we find four (Fig. H), three (Fig. |10|), two (Fig. [IlD, and finally one 'stripe' 
(Fig. [12]), where the infinite cluster spans the whole system. The spatial shape of stripes 
depends on lattice symmetry and boundary conditions. In the case of a two-dimensional 
square or triangular lattice with periodic boundary conditions, the system self-organizes 
into stripes with the boundaries along one of the principal axes. For absorbing boundary 



conditions, the stripes are replaced by regions of a different shape (see Fig. pT0|) . 

To understand the occurrence of stripes with different densities, we consider first a system 
with a very high density p <• 1. Since we are far above the percolation threshold p c ~ 0.593 
for random site percolation, the strength of the infinite cluster P is close to p. We start 
with a random initial state. In the first iteration step, we remove either the infinite cluster, 
which consists of nearly all occupied sites in the system, or one of the finite clusters, which, 
consequently, are very small. In the first case, only a few small clusters remain, and most 
of the occupied sites are redistributed randomly in the system. In the second case, only a 
small number of occupied sites are redistributed. In both cases, however, the state of the 
system changes little, and the new state is close to a completely random state. 

If we now decrease p, the remaining clusters (after the removal of the infinite cluster) 
become larger, and the density fluctuations increase. When the removed sites are refilled 
into the system, part of them are positioned in or near these surviving clusters, where the 
density then will be larger than the mean density. In the space between these clusters 
the density consequently is lower than the mean density (see Fig. O). If p falls below 



a certain threshold p*, this density between the surviving clusters becomes smaller than 
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the percolation threshold p c . Then, there exists no infinite cluster in the system after the 
first time step (or after a few iterations), and the state with one stripe becomes unstable. 
Evidently, p* > p c . 

The approximate value of p* can be derived by the following argument: The number 
of empty sites after the removal of the infinite cluster is L 2 (l — (p — P)). Let this subset 
of the total lattice be denoted by T>. It contains the area of the infinite cluster and the 
sites that have already been empty before. In the next step, the L 2 P occupied sites of the 
infinite cluster are redistributed randomly among these empty sites, leading to a density 
p' = P/(l — (p — P)) < p in T>. In the vicinity of the clusters that were left over the local 
density is now higher than before. One can think of this state approximately as consisting 



of many compact clusters embedded into a lower-density background (see also Fig. [L3|). It is 
then obvious that the homogeneous state with an infinite cluster stretching over the whole 
system cannot survive if p' < p c , since, in this case, the infinite network which links the 
finite high-density regions is broken. On the other hand, if p' > p c , the situation is not 
fundamentally different from the case p ^ 1. Fluctuations in the local density do exist, but 
they cannot get stronger with time, since the high-density regions are themselves removed 
from the system with high probability during the next few iteration steps. Thus, we arrive 
at the following implicit approximate equation for a threshold density p*. 

P{P = ^ = (10) 

1-P* 1-Pc 

In the simulations, we found p* ~ 0.625 for the square and 0.533 for the triangular lattice. 
A measurement of the left hand site of Eq. ([H]) at p = p* and L = 1024 yielded ~ 1.46 
for the square and m 1.01 for the triangular lattice. The value of the known right hand 
site is ~ 1.46 and 1, respectively. Although Eq. ([T0| ) is approximate, the agreement with 
simulation results is very good. 

If p falls below p*, the homogeneous phase becomes unstable, and the system rearranges 
itself to a new stationary state that consists of two homogeneous stripes with equal area and 
different densities pi and pi (pi > p-i). For p < p*, we have pi > p* and p2 < p* ■ When the 
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infinite cluster in this new state is removed, the density of stripe 1 decreases to p 3 = pi — P. 
Now these sites are re-injected randomly into the system, such that stripe 1 is filled up to 
P2 and stripe 2 to p\. The net effect is that the two stripes have changed their roles. 

If we lower p further, we will eventually reach the density p for which p\ = p*. At 
this point the two-stripe state becomes itself unstable (because the highest density stripe 
becomes unstable) and reorganizes into a three-stripe state, thereby increasing pi such that 
pi > p* again. The dynamics of the three-stripe state is analogous to the two-stripe state: 
When the infinite cluster is redistributed, the stripes simply exchange their densities (3 — > 2, 
2 — > 1, 1 — >• 3). With decreasing density, the number of stripes n increases further. 

The stripe structure breaks down when the system size becomes so small that the width 
of a stripe is of the same order as the roughness of its boundary. The resulting state shows 
patches of different size, but differs from the critical state in its cluster size distribution. 
In Fig. [L^, the cluster size distribution for a system at p = 0.46 in the 'patched' state 
(L = 512, 1024), and 'striped' state (L = 2048,4096) is plotted. The transition between the 
two seems continuous, due to finite-size effects. With increasing L, the bump becomes more 
distinct because the system tries to separate one infinite cluster from the ensemble of all 
clusters. There is no scaling of n(s), not even in the 'patched' states. With increasing system 
size, the bump should separate completely from the size distribution of finite clusters and 
move towards s = oo. The distribution of densities in the system changes also continuously 
with increasing L, when the transition from patches to stripes is made. Fig. [16] shows a 
histogram of the local density for L = 1024, i.e. for the patched state, averaged over 20 2 
sites. This histogram shows pronounced peaks that are precursors of the stripes. In the 
limit of infinite system size, these peaks will become infinitely sharp. 

The densities of the different stripes are related by several equations. Let p 1; . . . ,p n be 
the densities in a state with n stripes, starting with the highest density. Additionally, we 
define the density p n +i of the stripe which contained the infinite cluster, immediately after 
the infinite cluster has been removed from the system, and before the removed sites are 
refilled into the system. As consequence of a large redistribution, the different stripes just 
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exchange their densities, i.e. 

Pi-l = Pi + (pl ~ Pn+l) ■ (1 - Pi) I (n (1 - p) + pi - pn+l) (11) 

for i = 2, . . . ,n + l. The last factor on the right hand side represents the fraction of occupied 
sites of the infinite cluster that are refilled in the stripe with density pi. We finally obtain 

1 - Pi = 1 ~ P2 = = 1 ~ Pn ^ 
1 - p 2 1-/03 " 1 - Pn+l ~ 

Together with p = (1/n) Yn=iPi we have n equations for n + 1 densities. The average density 
in a system with n stripes is then 

l-p n+1 " / 1-pi v/ " 



, - Pn+l J- - PI 

p= — n — Mr^r 



Pn+l J 

1 Pl Z P n+1 (is) 

^(((i-po/a-Pn+i))- 1 /"-!)' 1 ) 



This equation was already used in Sec. IV 



The striped phase has several features in common with the critical and subcritical phase: 
They can all can be characterized by regions of different density. The density of these regions 
goes through a cycle, and the spanning cluster of a region is removed when its density is 
of the order p* w 62.5%. Although there is no strict synchronization in the subcritical and 
critical phase, the system seems to maintain many features of the synchronized phase, albeit 
with more irregularities and fluctuations. If there were no regions with pi oca i = P* \ they 
would be generated by the same mechanism as in the synchronized phase. The resemblance 
to the synchronized phase is especially evident in the one-dimensional case. There, the 
critical point can be interpreted as a synchronized state with p* = 1, p°° = 0, and an infinite 
number of stripes. 

From our understanding of the striped phase and its stability, we cannot rule out the 
existence of gaps, i.e. values of the density for which no synchronized phase is stable. The 
following scenario could occur: When decreasing the density of an n-stripe state below the 
stability threshold, the new density in stripe 2 after the restructuring could be too large (i.e. 
> p c ) for the state to be stable. In this case, the new {n + l)-stripe state would be stable 



only after lowering the overall density a bit further until P2 < p c - The intermediate state 
would then have some irregular patched structure. However, in our simulations we could 
not observe this phenomenon. 

B. Hysteresis and Maximum Number of Stripes 

Since the transition between states with different number of stripes is discontinuous, 
hysteresis effects are to be expected. When the overall density p is decreased adiabatically, 
a state becomes unstable when p x falls below p* ~ 0.625, and its number of stripes increases 
by one (thereby also increasing pi). When p is increased adiabatically, however, there is 
no reason, why the system should necessarily rearrange itself at p± — p*. The stripe with 
the highest density can never become unstable, since it only comes closer to the ideal case 
of -P(pi) = pi, where it consists only of the infinite cluster. The reorganization from an 
(n + l)-stripe state to an n-stripe state is then triggered by the stripe with second highest 
density p2- In the limit of infinite system size, it takes place when this stripe starts to 
contain an infinite cluster, i.e. when p 2 approaches the percolation threshold. As long as 
P2 < Pc, even the largest cluster in stripe 2 is small compared to the infinite extent of one 
stripe and cannot have any effect on the dynamics. For finite system size, however, the 
reorganization takes place as soon as the ratio $,2/ D exceeds a certain critical value, where 
D is a measure for the linear extent of one synchronized region (as e.g. the thickness of 
a stripe in the two-dimensional system with periodic boundary conditions). Therefore, the 
effect of hysteresis in a finite system is difficult to observe when the number of stripes n 
increases. Only for the transition from the two-stripe state to the homogeneous state we 
could identify a significant interval of hysteresis. The two-stripe state could be kept alive 
up to p = 0.64 for L = 4096. 

In the following, we argue that there exists an upper limit to the number of stripes n 
even in an infinitely large system. When lowering p, the difference p\ — P2 decreases with 
each additional stripe. If the number of stripes n was not bounded from above, we would 
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have pi — P2 — > and pi — > p*. Such a state could only be stable if the maximum cluster size 
was finite in all stripes but stripe 1, which, in turn, could only be possible if the percolation 
threshold was identical to p*. Since the distribution of occupied sites within a stripe is 
not completely random, the percolation threshold can be different from the threshold for 
conventional site percolation. We can find its value by simulating the density sequence 
(p n , p n -i, pi) for a single stripe. We start with a stationary state with density p just 
above p*. We remove the infinite cluster and fill the system again randomly and continually, 
and we measure the percentage of the largest cluster for two different system sizes. The 
results are shown in Tab. |J. We see from Tab. |I| that an infinite cluster exists already for 
p = 0.60, since P is non-vanishing and does not decrease with increasing L. Furthermore, 
we measured the correlation length of the finite clusters in the interval 0.57 < p < 0.62, 
which shows a pronounced peak around p = 0.59 for both L = 512 and L = 1024 (see 
Tab. U). These results show that the percolation threshold for stripes is very close (if not 
identical) to p c for site percolation. Thus, states with p 2 > 0.59 cannot exist, and there 
must exist a minimum density gap Ap = p\ — p 2 , and, consequently, a finite maximum 
number N of stripes. Eq. ( |T2"D gives a maximum possible number of n = 11 (±2) stripes, and 
a corresponding minimum mean density p = 42(±0.3)%. In our simulations, due to finite 
system size, we could observe a maximum number of n = 5 stripes at p = 0.45 and L = 4096 
(see Fig. [8]). For other lattice types, the maximum number of stripes and the densities where 
the phase transitions take place are of course different. The realization of an infinite number 
of stripes (i.e. a front moving through a continuum) is possible with different rules, where 
the occupied sites of the chosen cluster are removed one by one and are put back into the 
system, before the next site is removed (see [fl8[] ) . 



A schematic phase diagram of the system for all densities from to 1 is shown in Fig. [TJ. 
Since a state with n stripes is only stable if the thickness of a stripe is much larger than 
£2, the phase boundaries depend on the system size L for small L. They become vertical 
for large L, and the values of p, where the phase transitions take place, are well defined. 



As one can also see in Fig. |1J, the minimum density of the synchronized phase is smaller 
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than the upper limit p^f> of the critical phase. Since the transition from patches to stripes 
is discontinuous, it shows also hysteresis. 

C. Stability, Movement, and Roughness of Stripes 

In systems with n > 3 stripes, there exists the possibility that the stripes might not be 
arranged in consecutive order with respect to their densities. However, these arrangements 
are not stable, as will be shown in the following. If the stripes are arranged in consecutive 
order, the stripe with highest density has always the same densities in its neighboring stripes 
during the complete cycle of length n. Each time the infinite cluster is removed, the borders 
of the stripe it has occupied expand £ 2 to one side and £ n to the other. Since during one cycle 
the infinite cluster occupies each stripe position exactly once, all changes of width cancel. 
Since £ 2 > £n, the whole pattern moves into the direction of the stripe with second highest 
density, as seen from the infinite cluster. If the stripes are not arranged in consecutive order, 
e.g. {pi, p 3 , p 2 , Pi\, the changes in width do not cancel. The left neighbors of the infinite 
cluster during one cycle are p 4 , p 3 , p 3 , p 2 , and its right neighbors p 3 , p 2 , p 4 , p 3 . Therefore, there 
is a net change of the width of some stripes on the expense of others, and sooner or later 
(depending on system size) the structure breaks down, and the system rearranges itself to a 
state with the stripes arranged in consecutive order. The instability of certain configurations 
like the one mentioned above was tested by starting with an artificially generated stationary 
state with the 'wrong' order of stripes. After some time, the restructuring to the stable 
stationary state could be observed. 

Similarly, if the area of different stripes is different or if the values of the stripe densities 
are not the same after each rearrangement of occupied sites, different stripes see a different 
environment when they are removed. Such states can therefore no be stable. However, when 
they are close enough to the stable state, they return to it, as follows from the observed 
stability of striped states. In order to check explicitly the stability with respect to variations 
in the densities in a system with two stripes, we start with an unperturbed state at density 
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p° with two stripes of density p5 and p°, respectively. The strength of the infinite cluster is 
P° = P(Pi), and Pi — P° = pjj. Now we perturb the system by changing the densities to 
Pi — Pi + $ an d Pi = p\ — 8 with small 5 > 0. The overall density p° is kept constant. The 
strength of the infinite cluster is now P = P(pi) = P{Pi + 5) = P° + 5 ■ P'(p±). After the 
removal of the infinite cluster the densities are p% = p\ — P = p° — 5 ■ (P' — 1) and p 2 = p° — 5. 
Using Eq. ( |TT1) one obtains for the highest density in the next iteration step 

new _ _ c , / pO , c . pt\ 1 - P2 + j 

If this density is lower than the perturbed density pi = p? + 8, the system wants to return to 
the unperturbed configuration and the initial state is stable. For stability, thus, the following 
inequality has to hold: 

P2 d + {F+d P) 2-(pO + pO) + 5.P' <P2 + ^2-(pO + pO) +d - 

On the right hand side, we have used p? = p° + -P°(l — /0°)/(2 — (p° + p°)) from Eq. (^Tj). 
This inequality can be rewritten as 

P < 4 ' 1 -"°) +3P ° (14) 

(l_po /(2 (i_ p o + po )))(1 _ p o ) I > 



We measured all quantities that appear in Eq. (|14D in the density interval where the two- 
stripe state could be observed and found that it is completely contained in the interval where 
this condition is fulfilled. 

As already mentioned above, the whole pattern of stripes moves into the direction of 
the stripe with second highest density seen from the infinite cluster. This movement can be 
observed in the simulations for all states with n > 2. The net velocity in a two-stripe state 
is zero due to symmetry. We measured the velocity of this movement for stationary states 
with three and four stripes at p = 50% and 46%. The results are v ~ 4.3 ± 0.5 lattice sites 
per iteration step for p = 50% and v ~ 5.2 ± 0.5 for p = 46%. Since the degree to which the 
infinite cluster reaches into the neighboring stripes depends on their correlation lengths, the 



velocity is essentially determined by the density in stripe 2. As one can see from Tab. |TI 
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p2 is larger for p = 46% than for p = 50%, therefore the velocity is higher for p = 46%. The 
movement of the stripe boundaries can be seen in Fig. [T7|, where a sequence of stationary 
states of a three-stripe state is shown. 

As one can see in Fig. || the stripes seem to have rough rather than smooth boundaries. 
When the infinite cluster is removed, all finite clusters in the adjacent stripes that are 
connected to it are also removed. This leads to not only to a net velocity in one preferred 
direction, but also to a roughening of the stripe boundaries, even in situations where the 
(artificially generated) initial state consists of strictly horizontal stripes. 

The roughness of the interface is characterized by an exponent a that describes how the 
saturation width of the interface w sat scales with its length L, 

w sat (L) oc L a . 

The interface width w is defined as the root mean square of the differences y(x) — y(x) where 
y(x) (y(x)) is the (averaged) ^-coordinate of the interface (assuming horizontal stripes) 



w(L,t) 



with y{t) = (1/ L)J2i=iy(x,t). We assume that y{t) is single-valued, i.e. there are no 
overhangs. This assumption is always correct at sufficiently large scales. 

To characterize the time-dependent dynamics of the roughening of the interface width 
w(t), the growth exponent (3 is introduced by 

w(L,t) oc t 13 . 

For small t, there is no dependence on L, since the information on system size needs a certain 
time to spread. For an introduction to interface roughening, see [[19[ . 



We measured the roughness exponent a as well as the growth exponent (3 for stationary 
states with two and three stripes at p = 55% and 50%. States with a larger number of 
stripes can only survive in large systems, so we were not able to scan a broad enough range 
of system sizes L to determine an exponent. Since at sufficiently large scales there is an 
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up-down symmetry in the case n = 2, which is absent in the case n = 3, where the boundary 
moves in one preferred direction, we expect different exponents. 

The results for a and (3 are shown in Tab. [TIT]. Fig. ^ shows the roughening of the 
boundary for the cases n = 2 and n = 3. The simulation results are compatible with the 
values a = 1/2, (3 = 1/3 for n = 3 and a = 1/2, (3 = 1/4 for n = 2, although other 
values cannot be ruled out. The former exponents are those of the Kardar-Parisi- Zhang 
universality class [pQl , whereas the latter belong to the Edwards- Wilkinson universality 



class |[21| , describing an interface in thermal equilibrium. For all n > 3, the boundaries 



exhibit a net movement in one preferred direction, thus breaking the up-down symmetry. 
Therefore we expect to find the KPZ universality class also in the cases n > 4 that could 
not be investigated in the simulations. 

In the limit of a large number of stripes n — > oo, the movement of the infinite cluster is 
reminiscent of the movement of a front. To study the dynamics of such a front, we introduce 
a third state ('excited') that marks the sites belonging to the front. The excitation spreads 
at each time step to all of its nearest occupied neighbors leaving behind an empty site. 
In order to assure a constant density, one has to re-occupy randomly chosen empty sites. 
Starting with a flat excitation front in an initial state of average density p, after short time 
the front will either have disappeared or the system will have evolved to a stationary state, 
where the front propagates quasi-deterministically in one direction. After leaving the system 
at one end, it reenters at the opposite end. Although the front state can be interpreted as a 
stripe state with n = oo, there is one important difference. In the stripe state, the infinite 
cluster sees a medium with density below p c in front of it, whereas in the front state, this 
density is above p c . One can therefore not expect to find KPZ behavior as in the stripe 
state. Instead, the shortest paths (using only occupied sites) from the sites of the front at 
time t — to the sites of the front at some later time have a Gaussian distribution, leading 



to a width of the interface that scales as log(L), i.e., a = 0. Fig. [19] shows the saturation 
width as function of the system size, together with a logarithmic fit. 

One can also compare this behavior to the behavior of an excitation front in a percola- 
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tion system without periodic boundary conditions perpendicular to the front, i.e. a system 
of size L parallel to the front and of size 'oo' perpendicular to it. The 'infinite' extension 
perpendicular to the front would in practice be realized by providing a constant, homoge- 
neous density ahead of the front, such that the front will never pass through a region where 
it already has been before, and additional correlations due to sites left over from the last 



passing are eliminated. In this case, the stability arguments of Subsec. [V A| and of [|18 
which lead to the instability for p beforc < p* f are no longer valid, and one can investigate 
the properties of the front in the whole density interval of p beforc £ [p c ; ]_]. Then, one can 
also measure the exponent that describes how the velocity of the front vanishes when the 
density p before approaches p c . As in the previous paragraph, we expect that such a front has 
a roughness exponent a = 0. Models which yield a front of finite roughness must therefore 
be more complicated and include a dependence of the front propagation on the shape of 
the front, or memory effects. Such models, as the Kuramoto-Shivashinsky equation, or the 
model in [E2|, give a KPZ roughness exponent. 



VI. SUMMARY AND CONCLUSIONS 

To conclude, we have described a nonequilibrium percolation model which shows several 
new phenomena which are unknown in equilibrium percolation. Clusters of occupied sites 
are removed and refilled into the system at randomly chosen empty sites. For densities 
smaller than a critical density, there are only finite clusters, as in percolation theory. At 
the critical point, the critical exponents assume values different from percolation theory, 
and they do not satisfy a hyperscaling relation. For densities between the critical density 
~ 40.8% and a second density ~ 43.5%, the system remains critical, with the correlation 
length diverging slower than the system size. The value of the exponent that relates the 
correlation length with the system size depends on the density. For densities larger than 
~ 43.5%, the system has a finite number of regions of different densities. The number of 
regions depends on the density, and the transitions between states with different numbers 
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of regions are discontinuous and show hysteresis. The shape of the regions is striped for 
periodic boundary conditions. 

We were able to describe the dynamics of the system in the critical state as well as in the 
striped state in terms of patches of different densities. On length scales smaller than a 'finger 
thickness' (critical region) or the stripe thickness (synchronized phase), the density of a patch 
is fairly homogeneous and goes through a temporal cycle: Starting from a small density, the 
density increases until most occupied sites of the patch are removed at p ~ 62.5%, and the 
cycle restarts. The system is therefore synchronized on small length scales. In the critical 
interval, the density is not large enough to allow a synchronization over the total width of 
the system, leading to the observed power laws. A similar relation between critical behavior 
and incomplete synchronization was found in the earthquake model ||. There, complete 
synchronization is hindered by the boundary conditions. In our model, the global density 
conservation prevents synchronization over distances larger than some correlation length. 
In both models, the correlation length diverges slower than the system size. A model with 
similar phenomena (i.e., subcritical, critical, and supercritical phases) has been reported in 

El- 

The density cycles observed in our model cannot exist in equilibrium systems, which are 
invariant under time reversal. Therefore, fluctuations in an equilibrium system are usually 
small and well described by a Gauss distribution around some mean value (except for the 
neighborhood of a critical point where the next-order terms have to be added to the free- 
energy functional). One consequence is the extensivity of equilibrium systems, i.e. a part of 
an equilibrium system behaves like a smaller system of the size of the part. This is not true 
for our model. 

Although our simulations were performed in two dimensions, we expect the same general 
picture in higher dimensions: Unfortunately, in dimensions d > 2, one has to cope with severe 
finite-size effects. For example, the three-dimensional analogue of a two-stripe system would 
be a two-layer system. For this state to be stable, the correlation lengths in the layers have 
to be small compared to their thickness. With a given maximum number of sites ~ 1.6 x 10 7 
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that could be simulated, the maximum thickness of one of two layers in a three-dimensional 
system is only « 125 compared to ~ 2000 in two dimensions. Another problem is that 
a visualization of the states like in two dimensions is not possible. One has therefore to 
measure local densities in order to determine the number of stripes that coexist in a given 
state. For the three-dimensional system, we could verify the transition to two stripes at 



p* D »i 0.34. In Fig. |20|, we plotted a histogram of the local density in a system at p < pg D , 
proving the existence of two stripes. 

We also simulated a variant where not an arbitrary site is selected, but always a site 
that belongs to the largest cluster in the system. The simulations were carried out for 
system sizes up to L = 1024. Although the snapshots look different, since small clusters are 
no longer removed, we did not find different behavior. The important quantities, i.e., the 
critical exponents, p^\ p^ and p* remain the same. 

An interesting modification of the model does not conserve the density of occupied sites 
exactly, but only on the average. While such a modification would not change the behavior of 
an equilibrium system, we expect some important changes in our model. When the density 
is not strictly conserved, large global density fluctuations can occur, and a sharp distinction 
between a subcritical, a critical, and a supercritical phase is no longer possible. This new 
model will be treated in another publication |24] . 

The existence of these fundamental differences between equilibrium and many nonequi- 
librium systems makes it difficult to apply methods developed for the study of equilibrium 
critical behavior in nonequilibrium systems. It remains still a challenge to find a field- 
theoretical formalism that allows to analyze analytically the critical behavior of this and 
related models. 

S.C. was supported by the Deutsche Forschungsgemeinschaft (DFG) Contract No Schw 
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FIGURES 




FIG. 1. Stationary state at p = 39.3% and L = 1024 (periodic boundary conditions, square 
lattice with nearest neighbors, occupied sites are black, empty sites are white) 
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FIG. 3. Upper curve: cluster size distribution n(s) of a homogeneous system at the percolation 
threshold p c ~ 59.3%. Lower curve: cluster size distribution n(s) of a system with densities between 
p after and / o before ; distributed according to Eq. (||). Both systems have size L = 2048. The dashed 
straight lines have slope 2.06 and 2.15, respectively. 
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FIG. 4. Stationary state at the critical density p c ss 40.8% and L = 4096 (periodic boundary 
conditions, square lattice with nearest neighbors, occupied sites are black, empty sites are white) 
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FIG. 5. Stationary state at p = 43% and L = 4096 (periodic boundary conditions, square 
lattice with nearest neighbors, occupied sites are black, empty sites are white) 



34 



10 



10" 



10" 



10 



10" 





(a) 


V (b) 






i i i 



10 



10 10 

S 



10 10 

S/S 



MAX 



FIG. 6. Normalized size distribution of clusters for p = 0.43 and L = 512, 1024, 2048, 4096, 
(a) before and (b) after rescaling. 
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FIG. 8. A stationary state with five stripes at p = 0.45 and L = 4096 (periodic boundary 
conditions, square lattice with nearest neighbors, occupied sites are black, empty sites are white) 
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FIG. 9. Stationary state with four stripes at p = 0.46 and L = 4096 (periodic boundary- 
conditions, square lattice with nearest neighbors, occupied sites are black, empty sites are white) 
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FIG. 10. Stationary state with three regions of different density at p = 0.50 and L = 2048 
(absorbing boundary conditions, square lattice with nearest neighbors, occupied sites are black, 
empty sites are white) 
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FIG. 11. Stationary state with two stripes at p = 0.55 and L = 1024 (absorbing boundary- 
conditions, square lattice with nearest neighbors, occupied sites are black, empty sites are white) 
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FIG. 12. Stationary state at p = 0.63 and L = 1024 (absorbing boundary conditions, square 
lattice with nearest neighbors, occupied sites are black, empty sites are white) 
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FIG. 13. Schematic picture of the stationary state at a density just above p* . Light (dark) grey 
regions have lower (higher) than average density. 
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FIG. 15. Cluster size distribution of a system at p = 0.46 for various system sizes 
(L = 512,1024,2048,4096). 
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FIG. 16. Histogram of the local density piocal f° r L = 1024 and p = 0.46, averaged over 20 2 
lattice sites. 




FIG. 17. A sequence of stationary states of a three-stripe state at p = 50% and L = 1024. The 
snapshots were taken at times to, to + 51, to + 102, to + 153 and to + 204. 
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FIG. 20. Histogram of the local density piocal in a three-dimensional system for L = 100 and 
p « 0.34 in arbitrary units, averaged over ll 3 lattice sites. 
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TABLES 



p 0.41 


0.42 




0.43 


0.435 


01 1.46(8) 


1.72(6) 




1.86(8) 


1.92(8) 


4> 2 0.79(2) 


0.92(3) 




0.97(3) 


0.99(2) 


03 1-23(3) 


1.50(3) 




1.62(3) 


1.69(5) 


TABLE I. The exponents < 


^1,2,3 °f s max) an d S for various densities p^ 


< p < p { c ] (L = 512, 


1024, 2048, 4096) 










p 0.57 


0.58 


0.59 


0.60 


0.61 0.62 


P(L = 512) 






0.33 


0.46 0.52 


P(L = 1024) 






0.34 


0.47 0.53 


£(L = 512) 57 


85 


133 


62 


36 13 


£(L = 1024) 58 


89 


145 


75 


28 14 


TABLE II. The percentage P of the largest cluster and the correlation length £ for different 


densities p = 0.57, . . . , 0.62 with L = 512 and L = 


1024 






P 


0.46 




0.50 


0.55 


Nr. of stripes 


4 




3 


2 


P 1 


0.63 




0.66 


0.66 


P 2 


0.53 




0.51 


0.44 


a 






0.44(10) 


0.46(10) 


(3 






0.35(5) 


0.25(5) 


V 


5.2(5) 




4.3(5) 





TABLE III. Roughness exponent a, growth exponent (3, and velocity v of the stripes for various 



densities (200 < L < 2000). 
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